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Abstract 

We develop the formalism discussed previously in hep-ph/0601209 and 
hep-ph/0605246 to construct a kinetic theory that provides insight into 
the earliest "Glasma" stage of a high energy heavy ion collision. Particles 
produced from the decay of classical fields in the Glasma obey a Boltz- 
mann equation whose novel features include an inhomogeneous source 
term and new contributions to the collision term. We discuss the power 
counting associated with the different terms in the Boltzmann equation 
and outline the transition from the field dominated regime to the particle 
dominated regime in high energy heavy ion collisions. 
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1 Introduction 

In two previous papers [1,2], we introduced a formalism to compute multi- 
particle production in field theories coupled to strong time-dependent external 
sources. The QCD example of such a field theory is the Color Glass Condensate 
(CGC) [3-16]. For simplicity, we considered a cf> 3 theory; we believe however 
that most of our results are of general validity and can be extended to gauge 
theories [17]. 
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In this paper, we will address a problem in multi-particle production that was 
not considered in Refs. [1,2]. Specifically, the approach developed there did not 
include scattering processes that are important for the dynamics of the system 
at late times. These are the so called secular terms which are of higher order 
in the coupling constant (loop corrections) and are accompanied by growing 
powers of time [18-20]. The secular contributions must be resummed to obtain 
sensible results. In a quantum field theory, this resummation is performed in 
principle by solving the Dyson-Schwinger equations. In practice, the Dyson- 
Schwinger equations are difficult to solve. For a system of fields coupled to 
an ensemble of particles, it is well known that the Dyson-Schwinger equations 
can be approximated by a Boltzmann equation for the distribution of particles. 
The goal of the present paper is to extend the approach of Refs. [1,2] to derive 
a kinetic equation that includes the late time contributions to multi-particle 
production in field theories with strong external sources. We have in mind the 
dynamics after a heavy ion collision, where the classical field produced by the 
colliding nuclei expands rapidly into the vacuum along the beam direction. Our 
approach may also be of relevance to descriptions of the decay of the inflaton 
field and thermalization in the preheating and reheating phases of the early 
univcrsc-a nice review with relevant references can be found in Ref. [21]. In 
both cases, as the classical field evolves, the occupation number decreases and 
it is more appropriate to describe the higher momentum modes of the system 
in terms of particle degrees of freedom. 

The connections between the classical approximation in field theory and ki- 
netic equations in the framework of nuclear collisions were previously discussed 
by Mueller and Son [22], and subsequently by Jeon [23]. They considered a 
system of fields in the presence of an ensemble of particles described by a distri- 
bution /. Performing a classical approximation in the path integral describing 
the evolution of this system and a gradient expansion in the obtained Dyson- 
Schwinger equations, these authors obtained a kinetic equation for /. An obvi- 
ous question arises: with what accuracy does this kinetic equation reproduce the 
Boltzmann equation one would obtain without performing the classical approx- 
imation? The authors of Refs. [22,23] find that the kinetic equation obtained 
from the classical path integral reproduces correctly the collision term in the 
Boltzmann equation to leading power of / and (surprisingly) the first subleading 
term in / as well. 

We shall adopt a more ab initio approach here by considering a system that 
does not contain any particle degrees of freedom initially, but where the fields 
are coupled to a strong time-dependent external source j. The external source 
is assumed to be a stochastic variable that belongs to an ensemble of charges 
specified by a distribution W\j\. This is the typical set up in the description of 
heavy ion collisions in the Color Glass Condensate framework where W[j] repre- 
sents the distribution of color charges. Because of the expansion of the system, 
one may anticipate that the system can be described by field theory methods at 
early times and by kinetic theory and hydrodynamics at later times. The matter 
in this regime in heavy ion collisions has interesting properties; two noteworthy 
possibilities are dynamically generated topological charge [24,25] and plasma 
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instabilities possibly leading to turbulent color fields [26] . This matter has been 
called an Glasma [27,28] and understanding its dynamical evolution holds the 
key to a deeper understanding of the strongly interacting Quark Gluon Plasma 
(sQGP) that may be formed at later times [29]. The 

We will address here general questions about the dynamical evolution of such 
matter in the simplest possible context of a scalar (<f) 3 ) field theory 1 : 

i. What is the kinetic equation one obtains in Geld theories coupled to strong 
external sources? Knowing the answer to this question is important for one 
to handle correctly the transition region between a field theory description 
and kinetic theory. Indeed, one expects from the work in Refs. [22,23] that 
there exists a window in time where both approaches correctly describe 
the dynamics 2 . This suggests that the kinetic equation in the overlap 
regime must know about the coupling of sources to fields at earlier times. 
How is this manifest, how important is this effect and how does it go away 
? 

ii. What terms in the kinetic equation are important at different stages of 
the expansion? The previous question hints that we will obtain a kinetic 
equation that has additional terms absent in the conventional Boltzmann 
equation. We would like to understand how this generalized Boltzmann 
equation converges to the usual one at late times. 

The paper is organized as follows. In section 2, we shall remind the reader 
of relevant formulae in the derivation [1] of the average number (n) of produced 
particles. In section 3, we shall write down the Dyson-Schwinger equations for 
the two-point functions in theories with time dependent strong sources. These 
provide the starting point for a derivation in section 4 of the corresponding 
kinetic equation for the Glasma. We observe that the coupling of the field to 
an external source leads to an inhomogeneous term in this kinetic equation. 
In section 5, we discuss the properties of the different terms appearing in the 
kinetic equation. Albeit the collision term in the kinetic equation looks iden- 
tical to the collision term in the usual Boltzmann equation, it contains novel 
contributions to the self energy that are of 0-loop and 1-loop order. We discuss 
the power counting for these different contributions and assess their relative 
contribution at different stages of the temporal evolution of the Glasma. We 
conclude with a brief summary and outlook emphasizing unresolved issues. An 
appendix addresses how the averaging over the sources j in our formalism can 
be re-expressed in terms of the usual ensemble average implicit in the derivation 
of kinetic equations. 

1 Even the "simple" scalar theory is non-trivial. It will indeed contain very general fea- 
tures of relevance to the Glasma albeit the latter will have significant (and very interesting) 
additional features that are absent in the scalar case. 

2 This has to be the case if one wants the final result to be independent of the time at which 
one switches between the two descriptions. 
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2 Ab initio computation of (n) 

We consider the theory of a real scalar field <p with cubic self-interactions, cou- 
pled to an external time dependent source j(x). The Lagrangian of the model 
is 

C ee \d»m - V - + j<f> . (1) 

In [1], we systematically calculated particle production from these sources. In 
the Color Glass Condensate framework that this toy model mimics, the colliding 
projectiles are represented by a statistical ensemble of currents j. Physical 
quantities are obtained by averaging over all possible realizations of the j's. In 
this section, we shall discuss the calculation of the average number of produced 
particles in a given configuration of j's. 

A general formula for the average number (n) of produced particles is 

^ = / (2^*2E P <°in|«Lt(p)«out(p)|O in ) . (2) 

The number of particles produced with a certain momentum p is defined as 
the expectation value of the "out" number operator in the initial state. This 
formula gives the number of particles at asymptotic times, after the particles 
have decoupled 3 . 

A simple reduction formula gives [30] 

<0 in |aL(pKut(p)|0 in > = | J d^xd^y e^e?™ 

x(a x + m 2 )(a y + m 2 ) (O in |0(z)<Xt/)|Oi„) , (3) 

where Z is the wave function remormalization factor. The expectation value 
in the right hand side of this equation has two important features : (i) the 
vacuum state is the "in" vacuum state on both sides and, (ii) the two fields 
inside the correlator are not time-ordered. The Schwingcr-Keldysh formalism 
[31,32] provides techniques for computing these types of correlators. 

The operators □ + m 2 amputate the external legs of the two-point function 
G-+(x,y) ee (0in\<t>(x)<t>(y)\0in) ■ Defining 

G-+ (x, y) ee — y - G-+ (x, y) , (4) 

we can write the average multiplicity as 

d 3 p 



Introducing the variables 



(5) 



v x + y la\ 

X EE —— , TEEX-y , (6) 



3 The "number of particles" at some intermediate time, while the fields arc still interacting, 
is not a well defined concept. 
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we can rewrite this formula as 



E >tr = i^'" ,v<; ■■ iV "' 



where 



G-+ (X, p)^J d\ e-^ G. + [x+ r -,X- r -) 



(7) 

(8) 



is the Wigner transform of G ^(x, y). 

In the Schwinger-Keldysh formalism, the propagators G ee '{x,y), (e, e' — 
+ , — ) can be expressed as 



G ee ,(x,y) = — 



iSj e (x) i5j e ,(y) 



e iV SK [j+,j- 



(9) 



where iV SK \j+,j-] is the sum of all connected vacuum- vacuum diagrams. 
When j + = j_ = j, iV SK [j, j] = and the sum of all vacuum- vacuum dia- 
grams is unity. 

Working out the functional derivatives, 



G ee /(x,y) 



SiV SK 5iV SK 



S 2 iV a 



i5j e {x) i5j e '{y) i5j € (x)i5j € >(y) 



(10) 



j +=]-=] 



As iV SK is the sum of connected vacuum-vacuum diagrams, any of its deriva- 



x y , x y 
— • • — - + • — - 



Figure 1: Diagrammatic representation of the disconnected (left) and connected 
(right) terms in eq. (10). The gray blobs denote the remnants of Green's func- 
tions after the free propagators at the endpoints are amputated. 

tives with respect to j± is a connected Green's function. Therefore, G ee ' can be 
decomposed as 

G ee ,(x,y) = G c €€ ,(x,y) + G™,(x,y). (11) 
These are, respectively, the connected part 



i6j e (x)iSj e >(y) 



(12) 



j+=j- 



and a disconnected part corresponding to the product of the expectation values 
of the field at the points x and y : 



G n J,(x,y) = (4>{x))(<t>(y)) with (<f>(x)) 



i5j±(x) 



(13) 



J+=3- 
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Figure 2: Example of a tree diagram contributing to the field expectation value. 
The black dots terminating branches of the tree represent insertions of the source 
j in the diagram on the right. The sum of these tree diagrams is represented 
(left) by a line attached to a gray blob. 

When j + = j- = j, the expectation value of the field is the same on the upper 
and lower branches of the contour: ^ + (x)) = (<^>_(a;)). This explains why we 
omitted the +/— index in the expectation value of the field. 

A typical tree- level contribution to (<j>(x)) is shown in figure 2. Note also 
that {<p{x)) vanishes if the external source j(x) is zero 4 . At tree level, because 
j + = j- = j, the sum over the +/— indices in the Schwinger-Keldysh formalism 
at all the internal vertices of the tree (including the sources) can be performed 
by using the identities 

Gl + -G° + _ = Gl , G°_ + -G __ = G° R , (14) 

where G° R is the free retarded propagator 5 . When this sum is performed, all 
propagators in the tree diagram can be simply replaced by retarded propagators. 
This is equivalent to the statement that (4>(x)) is the retarded solution of the 
classical equation of motion, 

(n + m 2 )<P(x) + ^<p 2 (x) =j(x) , (15) 

with a vanishing boundary condition at x n = — oo. 

Eq. (7) is the complete answer to the problem of particle production in the 
effective theory described by the Lagrangian of eq. (1). If one were able to 
compute G(x, y) to all orders, this formula would contain everything one needs. 
There would be no need for tools such as kinetic theory. 

4 We assume that the self-interactions of the fields are such that there is no spontaneous 
breakdown of symmetry when j = 0. 

5 In momentum space, this propagator reads G R (p) = i/(p 2 — m? + ip°e). 
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However, evaluating eq. (7) to all orders is an unrealistic goal. What has 
been implemented thus far is the evaluation of eq. (7) at leading order (tree 
level) to calculate the gluon yield in high-energy nucleus- nucleus collisions [33- 
38]. In [1], an algorithm was sketched to compute at next-to-leading order 
(one loop) in terms of the retarded classical field and of retarded fluctuations 
propagating in the classical field background. 

In practice, one has to truncate the loop expansion. As we will discuss in 
the next section, the correct way to perform practical calculations is within the 
framework of the Dyson-Schwinger equations. 



3 Dyson-Schwinger equations 

The main problem with the loop expansion described in the previous section 
is that, in general, truncations in G(x,y) will lead to an incorrect large time 
limit of the number of produced particles. This can be traced to secular terms 
containing powers of the time that invalidate the perturbative series in the large 
time limit. This can be cured by appropriate resummation; the well known way 
to do this is to solve Dyson-Schwinger equations [18-20]. In this section, we 
shall discuss the Dyson-Schwinger equations obeyed by the two-point functions 
Gee 1 {x, y) of the Schwinger-Keldysh formalism. We will see that the presence 
of a disconnected contribution to these 2-point functions leads to interesting 
features in the corresponding Dyson-Schwinger equations. 



3.1 Dyson-Schwinger equation for the connected part 

It is straightforward to write a Dyson-Schwinger equation for the connected part 
of the 2-point function, G\ fJ , that resums self-energy corrections : 

G c (x,y) = G°(x,y) + J d A u d A v G°{x, u) [ - iS(u, v)] G c {v, y) , (16) 

where — iS is a 1-particle irreducible connected 6 self-energy, evaluated in the 
presence of external sources. We shall not write here explicitly the ± indices car- 
ried by the various objects. Instead, we write the time integrations as integrals 
over the complete Schwinger-Keldysh contour C. 

It is convenient to extract from this self-energy a local piece, by writing 

H(u,v)=g&{u)6(u-v) + Tl(u,v). (17) 

Except for the background field, which is a genuine local contribution to the 
self-energy, there is a certain arbitrariness in this separation because it depends 
on the momentum scale at which we resolve the system. A contribution to the 
self-energy that does not change significantly over space-time scales on the order 
of the Compton wavelength p^ 1 can be treated as a mean field at that scale. 



6 It is connected in order to have a connected 2-point function after the resummation and 
it needs to be 1PI to prevent double counting. 
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O = j + o 



Figure 3: Top: Diagrammatic representation of the Dyson-Schwinger equation 
of eq. (16). The large gray blob denotes the 1-particle irreducible 2-point func- 
tion S. Bottom: decomposition of the 1PI self-energy £ into a local part g& 
and a non-local part II (denoted by a large light-gray blob), following eq. (17). 

Therefore, the mean field term 3>(u) will contain the classical field, and possibly 
changes in the dispersion relation due to medium effects 7 . 
The Dyson-Schwinger equation then becomes 

G c (x,y) = G°(x,y)-ig J^d 4 u G (x,u)<f>(u)G c (u,y) 

d i ud 4 vG°{x,u) -iU.{u,v) G c (v,y) . (18) 

Using 

p x + m 2 ]G°(x,y) = -i5 c {x-y) , (19) 

where S c denotes the delta function on the closed time path 8 , we can rewrite 
this equation as 

[D x +m 2 + g&(x)] G c (x, y) = -id c {x - y) - j d 4 u U(x, u)G c (u, y) . (20) 

3.2 Dyson-Schwinger equation for the disconnected part 

We also need a Dyson-Schwinger equation for the disconnected part of the 
Green's function, 

GT(x,y) = (<f>(x))(<t>(y)) ■ (21) 

Because the expectation value (0) is a connected 1-point function, it is natural 
to factor the connected propagator out of it, by writing 

(<t>{x))= J d i uG c {x,u)S{u) , (22) 



7 To allow for this possibility, we denote the mean field piece by a symbol distinct from the 
one used for the classical field. 

8 <5 C (x — y) = unless x° and y° are equal and lie on the same branch of the time path. 
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where S(u) is an "effective source" term 9 . By construction, one obtains 

p x + to 2 + g&(x)] (<j)(x)) = -iS(x) - J d 4 u n(x, u)(0(«)) . (23) 
Multiplying both sides by (<^(j/)), one obtains 
[□ x + m 2 + ff *(a;)] G nc (x, y) = -iS(x)(<f>(y)) - J d A u U(x, u)G nc (u, y) . (24) 
Defining 

-m s (x,y) = S(x)S(y) , (25) 
we can rewrite this equation as 

p x + m 2 +g$(x)]G nc (x,y) = - J d 4 u [iF (x, u)G c (u, y)+U(x, u)G nc (u, y) . 

(26) 

Adding eqs. (20) and (26), we obtain the Dyson-Schwinger equation for the 
complete two-point function: 

[D x +m 2 + G(x, y) = -iS c (x - y) 

-J d 4 u n s (x,u)G c (u,y) + U(x,u)G(u,y) .(27) 

The only formal difference between this Dyson-Schwinger equation and the equa- 
tion one obtains in the absence of the source j is the term proportional to II s 
in the right hand side. 

In principle, the resummations performed by solving eqs. (20) and (26) (or, 
equivalently, eqs. (20) and (27)) would completely cure the problem of secular 
terms. Such an approach has been pursued numerically in [39], but has not been 
attempted yet in the context of heavy ion collisions in the CGC framework. 



4 Kinetic equation 

The Dyson-Schwinger equations we wrote down in the previous section contain 
all the necessary physics but their solution is likely too difficult; they there- 
fore by themselves do not provide any practical insight into the dynamics of 
high energy heavy ion collisions. One can simplify the problem a step further 
by transforming the Dyson-Schwinger equations for the 2-point functions into 
kinetic equations. However, as we shall discuss shortly, doing so requires that 
certain assumptions be satisfied. 

9 In the classical limit, one has S(x) = j(x) + £<j> 2 (x) (see section 5.2). 
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4.1 Fields and particles 

As is well known, the Boltzmann kinetic equation describes the space-time 
evolution of particle phase space densities. Therefore, to achieve a kinetic de- 
scription, the formalism considered thus far should be extended to incorporate 
an ensemble of particles. This is simply done by modifying the free propagators 
to add a term that depends on the distribution of particles f(p). In momentum 
space, the modified propagators are 10 

G° ++ (p) ee - L— - + 2nf(p)S(p 2 m 2 ) , 

p z - m l + ie 

G __(p) ee =i , + 2nf(p)S(p 2 - m 2 ) , 
p z — m A — it 

G°_ + (p) ee 2tt(%°) + f(p))S(p 2 - m 2 ) , 

G° + _(p) ee 2n(9(-p°) + f(p))5(p 2 - m 2 ) . (28) 

These modified rules for the Schwinger-Keldysh propagators can be de- 
rived [40] when the initial density matrix that describes the ensemble has the 
form 



exp 



/Jo 
{ 2n/ 2Ep PP E P ^n(P)«in(p)] , (29) 



where (3 P is a momentum dependent quantity. (Note: /3 P should not be confused 
with the inverse temperature.) Such a form for the density matrix is required 
if correlators computed with this density matrix are to satisfy Wick's theorem. 
From this form of the density matrix, one obtains the Schwinger-Keldysh rules 
of eqs. (28), with 

The function f(p) in the propagators only represents the initial distribution 
of particles in the system. Thus the field theory defined by the Lagrangian of 
eq. (1) and the propagators of eqs. (28) describes a system of fields coupled to 
an external source j and to an ensemble of particles with an initial distribution 
f(p). The Feynman rules then enable one to calculate the properties of this 
system at a later time. 

However, eqs. (28) do not lead to a well behaved perturbative expansion, 
except when the function f(p) is the equilibrium Bosc-Einstein distribution 
in our model of bosonic fields. In general, when f(p) is not a Bose-Einstein 
distribution, the perturbative expansion based on eqs. (28) is plagued by the 
previously mentioned pathological secular terms which need to be resummed. 
The time-scale at which resummation becomes necessary is related to the trans- 
port mean free path in the system, namely, the time between two large angle 
scatterings undergone by a particle. This resummation makes the distribution 
f(p) time-dependent reflecting the changes induced by collisions on the particle 



10 Thc propagators of the Schwinger-Keldysh formalism appropriate for calculating eq. (3) 
are the same with f(p) = 0. 
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phase space distribution. Under certain approximations to be discussed later, 
this temporal evolution is governed by a Boltzmann equation. 

The problem formulated in section 2 concerned a system that has no ensem- 
ble of particles at the initial time (f(p) — in eqs. (28)). At first sight, as / = 
is a particular case of the Bose-Einstein distribution (with a vanishing temper- 
ature), secular divergences may appear to be absent. However, this conclusion 
is incorrect because of the presence of external sources which drive the system 
out of equilibrium. Thus it is also necessary to resum secular terms in this case, 
leading to changes in f(p). The generalized propagators in cq. (28) constitute 
the natural framework to achieve this. Because the external source is both time 
and space dependent, one has more generally 

f(p) - f(X,p) (31) 

in eqs. (28). 

An important point must be made here about the tree level expectation val- 
ues (4>±(x)) in this f-dependent extension of our formalism. A crucial property 
of the propagators in eqs. (28) is that they still obey eqs. (14). The retarded 
propagator is therefore f -independent. Therefore, as long as loop corrections 
are not included, the field expectation value does not depend on / and is iden- 
tical to the result obtained from the retarded solution of the classical equations 
of motion. Hence, the contribution from the disconnected part of the 2-point 
function lead to an inhomogeneous (/-independent) term in the Boltzmann 
equation. 

4.2 Gradient expansion 

The extension (28) of the propagators leads to Dyson-Schwinger equations that 
are formally identical to eqs. (20) and (26) - with all the building blocks now 
constructed with /-dependent propagators. The first step in obtaining the 
Boltzmann equation is to rewrite all the distributions in terms of their Wigncr 
transforms. For a two-point function F(x,y), its Wigner transform F(X,p) is 
defined to be 

F(X,p) = J d 4 s e-*' F (x+ S -,X - |) . (32) 

The next step is to perform a gradient expansion where only long wavelength, 
low momentum modes are retained. In particular, all terms of order two or 
higher in d x are neglected. As our goal is to construct a kinetic theory for the 
Glasma, we will discuss the validity of this gradient expansion in the context of 
heavy ion collisions in the CGC framework. In this framework [14-16], the color 
sources p a (x±) generating the color currents 11 are stochastic variables that 
vary from event to event with a distribution H^p]. When calculating a given 
physical quantity, one first computes it for an arbitrary p and then averages over 
all possible p's in the ensemble generated with the weight W[p\. For example, 

11 These color sources are the QCD analogs of the sources j in our toy scalar theory. 
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in the McLerran-Venugopalan model [3-5], the distribution W[p] is a Gaussian 
with 

W[p] = exp ( - / * x± * v± , (33) 



where 



2 (x_L,y_L) 



p 2 (x ± ,y ± ) = (p(x ± )p(y ± )) = p 2 A (x ± )S(x ± - y ± ) . 



(34) 



Here p A (x±) represents the density of color charges at a spatial position x±_ 
in the nucleus. The typical momentum scale of the sources-the saturation mo- 
mentum squared Q 2 S at x± is simply related to p? A . 

The difference between one particular element of the ensemble and the aver- 
age weighted by W[p] is illustrated in figure 4 for the quadratic form p 2 . Because 





Figure 4: Left: p 2 distribution for one configuration in the ensemble represented 
by the distribution W[p\. Right: ensemble average of (p 2 )- 



the p(x±) are uncorrelated at different points in the transverse plane of the nu- 
cleus, a particular configuration of p's leads to a very rough density profile; in 
contrast, the average smoothly follows the Woods-Saxon density profile of a 
nucleus. This example simply illustrates that the gradients are uncontrollably 
large for a given configuration p rendering any gradient expansion meaningless. 
On the other hand, it is perfectly legitimate for ensemble averaged quantities. 

The typical momenta of "hard" particles is set by the saturation scale which 
is of order Q s ~ 1-2 GeV at RHIC energies; this scale may be higher at the 
LHC. In contrast, the gradient d x for averaged quantities changes appreciably 
over distance scales of the inverse nuclear radius given by ~ R^ 1 ~ 40 MeV for 
a large nucleus. The small magnitude of this scale in the gradient expansion 
relative to the typical saturation momentum justifies the gradient expansion for 
quantities that are averaged over the ensemble of color charges. 

The corresponding changes to the Feynman rules are described in appendix 
A. Here it is sufficient to note that the ensemble average is obtained by con- 
necting all the external sources j in the manner specified by the distribution 
VF[p]. For instance, in the MV model W[p] is a Gaussian, which implies that all 
the sources must be connected pairwise. The objects II s , II and <fr that appear 
in the Dyson-Schwinger equations (20) and (26) must be thought of as being 
averaged over j. In Feynman diagrams, we will represent the average over j by 
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surrounding the diagram by a light gray halo 



(35) 



This compact notation encompasses a very large number of contributions. For 
instance, at leading order, one would first approximate (^(x)) as the sum of all 
the tree diagrams, an example of which is represented in figure 2. For each such 
tree diagram, the sources j (the black dots in figure 2) are reconnected pairwise 
in all the possible ways. A typical reconnection of the sources, corresponding 
to the topology of figure 2, is displayed in figure 5. Note that the "loop order" 



Figure 5: Example of a tree level contribution to the average over the sources j 
of the field expectation value for a Gaussian distribution of sources. The links 
in red represent the elementary correlators (j(a;)j(j/)). The source connections 
represented here are for simplicity among nearest neighbors; all other pairwise 
topologies are feasible. 

of a given diagram is a meaningful concept only for diagrams before they are 
averaged over j. Indeed, as one can see by comparing the figures 2 and 5, the 
diagram before the j-average has loops and is of order g~ l . After the average 
is performed, while it has a large number of "loops" which do not contain any 
information about the order in g of the diagram. 

4.3 Boltzmann equation 

The final ingredient in the derivation of the Boltzmann equation is the so-called 
"quasi-particle ansatz" which can be expressed as 

G_ + (X,p) = (1 + f(X,p))p(X,p) , 

G+-(X,p) = f(X,p)p(X,p), (36) 
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where the spectral function p(X,p) is 

p(X,p) ee G R (X,p) - G A (X,p) = G-+{X,p) - G+-{X,p) . (37) 

The physical assumption here is that the interactions in the system are such 
that the collisional width of the dressed particles remains small compared to 
their energy; the system is made up of long-lived quasi-particles. 
The Boltzmann equation can now be obtained as follows: 

i. Write a Dyson-Schwinger equation analogous to eq. (27), but with the 
differential operator □ + m 2 + g& acting on the variable y instead of x, 
and subtract it from eq. (27). 

ii. Rewrite this equation in terms of the Wigner transformed quantities and 
perform a gradient expansion keeping only leading terms in d x . 

iii. Replace the Green's functions with the quasi-particlc ansatz and drop the 
spectral function p(X,p) which appears as a factor in all the terms. 

If the terms proportional to II s were absent from eq. (27), the steps outlined 
above would result in the well-known Boltzmann- Vlasov equation, 

2p ■ d x f(X,p)+gd x &(X) ■ d p f(X,p) = 

= (1 + f(X,p))U + -(X,p) - f(X,p)Tl- + (X,p) . (38) 

The extra term we have in the Dyson-Schwinger equations, proportional to II s , 
will modify the Boltzmann- Vlasov equation. Two key features of this novel term 
will prove essential in our derivation. The first is that II s (x, y) does not depend 
on whether the points x and y are on the upper or lower branch of the time 
contour. This is because the expectation value of the field, for equal values 
of the sources j+ and j_, is the same on both branches of the contour. The 
second feature is that the non-connected part of the propagators drops out of 
the spectral function, for the same reason. Hence, 

p(X,p) = G_+(X,p) - G+_(X,p) = G c _ + (X,p) - G c + _(X,p) . (39) 

Utilizing these two properties, we can perform the gradient expansion for this 
additional term in the same way as performed for the usual self-energy correc- 
tion. It modifies the right hand side of the Boltzmann equation by an additive 
correction 12 II s (X,p). Therefore, our final expression for the kinetic equation 
is 

2p-d x f(X,p) + gd x *(X) ■ d p f(X,p) = 

= U s (X lP ) + (1 + f(X,p))U + _(X,p) - f(X,p)Tl_ + (X,p) . 
(40) 

12 Notc that prior to dropping the spectral function that appears in all terms, we would have 

n s + _ (x, P )Gl + (x, P ) - n s _ + (x, P )G c + _(x, P ) = n« (x, P ) [gi + (x, p ) - g c + _(x, p )] 

= ns(x,p) P (x,p) . 
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The novel "source term" II s (A, p) in this equation is non-zero even if the particle 
distribution f(X,p) is zero. It is therefore responsible for / = not being a 
fixed point of the above equation; the solution of this equation is non-zero at 
later times even if the initial condition had a vanishing particle distribution. 
In the next section, we will discuss further significant differences between this 
kinetic equation and the conventional Boltzmann-Vlasov equation in eq. (38). 

5 Properties of the Glasma kinetic equation 

In this section, we shall discuss the various terms in eq. (40) with emphasis on the 
differences between these and those appearing in the conventional Boltzmann 
kinetic equation. 

5.1 Vlasov term 

We first consider the Vlasov term {gd x <f> ■ d p f) in the Boltzmann equation. We 
note that in performing the average of the mean Geld <I?(A), over the external 
sources j, the various correlation functions (j(xi) ■ ■ ■ j{x n )\ permitted by the 
distribution of sources W[j] are nearly translation invariant. The dependence 
of these correlators on the barycentric co-ordinate X = (x± + ■ ■ ■ + x n )/n is very 
slow because it arises from the density profile of the colliding nuclei 13 . Therefore 
the 1-point function 3?(A), averaged over j, also has a very slow dependence on 
its argument X; its Fourier transform with respect to X has only modes with 
momenta on the order of the inverse nuclear radius. As discussed previously, 
this scale is very small relative to the typical momentum of the particles under 
consideration and it is therefore legitimate to approximate it as a Vlasov term. 

As is well known, the effect of this term in the Boltzmann equation is to 
change the momentum of particles as they move between regions where the ex- 
ternal field is different. Indeed, gd x <& is the force that acts on the particles at 
point X and accelerates them towards regions of lower potential . The mean 
field <I> includes not only the classical field directly produced by the external 
sources, but also possibly a contribution coming from the particles encoded in 
f(X,p). Such a modification may arise from a modification of the particle dis- 
persion relation due to the collective action of the other particles. For instance, 
if the particles acquire a medium mass with a weak space-time dependence, 
this mass can be represented by a potential in the Vlasov term of the kinetic 
equation. 

13 The fact that this density profile is not a constant is the only effect in the problem that 
breaks translation invariance. 

14 For non central collisions, the shape of the overlap region between the two nuclei is elliptic; 
one has stronger gradients in the direction of the small axis of the ellipsis relative to those in 
the direction of its large axis. The Vlasov term therefore accelerates particles preferentially 
in the direction of the small axis of the overlap region. This leads eventually to elliptic flow 
and to an anisotropy of the spectrum of particles in momentum space. This effect is obtained 
entirely within kinetic theory without any assumption about the degree of thermalization of 
the system. 
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5.2 Source term in the kinetic equation 

Let us now consider the effect of the source term II s (X, p) in eq. (40), which can 
be obtained as the Wigner transform of the product S(x)S(y). An interesting 
situation, relevant for heavy ion collisions, is when tree diagrams are dominant 
because the external source is strong (gj ~ 1). In this case, the expectation 
value (4>{x)) is dominated by the retarded classical field (f>(x); the connected 
part of the 2-point function, G c , is simply the propagator of a fluctuation on 
top of the classical field, 

(GT 1 = D + m 2 + #. (41) 
One therefore immediately obtains the following expression 15 for S(x) : 

S(x) = □ + m 2 + gcj){x) 4>{x) 

= j(x) + 9 -4> 2 (x). (42) 

We see here that the effective source S(x) receives two contributions : 

i. the external source j(x) itself. This term is only important if we want to 
use the Boltzmann equation in regions of space-time where the external 
source is still active. In a heavy ion collision, the color sources are present 
only on the light-cone at a proper time r = 0. We will not consider this 
term further. 

ii. A term quadratic in the classical field produced by the external source; 
this term continues to contribute after the external sources have stopped 
acting. 

One may represent this effective source graphically as 

S(x)= -. = + . (43) 



The second term has a fairly straightforward interpretation. When the term 
quadratic in (f> in the classical equation of motion 



<j>{x)=j{x)- 9 -<t> 2 {x), (44) 



is important, we see that the field is not a free field. If expanded in particle 
modes, the number of particles in the field would change with time. Therefore, 
if one switches between a description in terms of classical fields to the kinetic 



15 This is the result for a potential g<f>^/3\. For an arbitrary potential V(<j>), the expression 
of S(x) in this approximation would read 

S(x) = j(x) - V'Wx)) + Hx)V"(4>(x)) , 
where the prime denotes a derivative of the potential with respect to tj>. 
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equation at a stage where this non-linear term is still significant, the source 
term in the Boltzmann equation modifies the number of particles in order to 
take this effect into account. 

At tree level, the effective source S(x), and hence II s , is independent of 
the distribution of particles /. As discussed previously, this is a straightforward 
consequence of the fact that, at tree level, the 1-point function in the Schwinger- 
Keldysh formalism can be rewritten entirely in terms of retarded propagators 
that are /-independent. II s is therefore non-zero even if / = 0. In contrast, 
the terms II± T in the r.h.s of the Boltzmann equation depend on / and vanish 
when / = as expected for collision terms. II s is therefore a source term in 
the Boltzmann equation, because it drives / to a non-zero value even if one has 
/ = initially. 

When we perform the average over j of the disconnected product S(x)S(y), 
we get both disconnected and connected source terms, 



depending on how the sources j are reconnected. In this picture, each light 
shaded area is simply connected after the average over j has been performed, 
and all the sources j it contains are linked in all the possible ways that preserve 
its connectedness. 

The first term in the r.h.s. of eq. (45) corresponds to contributions where 
we connect together only j's that belong to the same factor S, (S'(x))_ ) .(5(j/)) j .. 
Our previous remark about the average over j of the 1-point function &(X) also 
applies here to (S'(x))^ : its Fourier transform only contains very soft modes 
of the order of the inverse of the nuclear radius. It is therefore nearly zero for 
the typical particle momentum p ~ Q s we are interested in here. Thus only the 
connected terms in the average of the source term (S(x)S(y)) . matter in the 
kinetic equation. 

5.3 Magnitude of field insertions 

The source term in eq. (45), as well as the other terms in the right hand side 
of the Boltzmann equation, involve insertions of the classical field 4>(x). In this 
subsection, we present a simple power counting that enables us to estimate the 
magnitude of such insertions. To simplify the discussion, we shall assume that 
the space-time coordinate X corresponds to sufficiently late times when the 
external source j is zero and its influence is only felt through the classical field 
4>(X) generated by the source at earlier times. 

Following the discussion after eq. (34), we assume that there is hard mo- 
mentum scale Q s in the problcm-the saturation scale in heavy ion collisions. 
Typical particle momenta are of order p ~ Q s . In our toy model, the coupling 
constant g has the dimension of a mass in 4 dimensions. To mimic the power 




(45) 
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counting in QCD, we will write it as 



g = XQ s , (46) 

where A, like the QCD coupling constant, is dimensionless. We assume that the 
coupling constant Acl. 

To estimate the order of magnitude of the source term given in cq. (45), it 
is not sufficient to know the magnitude of the classical field. Kincmatical phase 
space constraints can alter the naive power counting. As these considerations 
will apply equally to the collision terms in the Boltzmann equation, it is worth 
our while to discuss the power counting for the source term at length here. 

From eq. (25) and eq. (45), the naive power counting for the source term 
would give 

IF(X,p) = ^y dhe^ ^ { X+ S -) < p 2 (X- S -)) j . (47) 

We will demonstrate that eq. (47) vanishes when the momentum carried by 
the classical field <p is nearly on shell. Rewriting this expression entirely in 
momentum space in terms of the Fourier transform <j>{k) of the classical field, 

IF(X ») = ^i [d 4 s f^L...^± e ^ s 
[A ' P) - 4 J dS J (2tt) 4 (2tt) 4 6 

x e - lfcl .(x+ f ) e - l fe 2 .(x+f) e - 4 fc 3 .(x-f) e - 4 fc 4 .(x-f) (^(kS(k 2 )0(kS(k 4 )) . 

(48) 

For the sake of simplicity, let us assume that the average over the external source 
j of the product of four fields factorizes into products of averages of two fields 
as suggested by the source distribution in eq. (33). 



(/>(fci)0(fc 2 )</H&3)<M&4)/ = (^(^1)0(^3)/ (^(^2)0(^4)/ +othcr contractions . 

(49) 

For illustrative purposes, we consider only one of the possible contractions corre- 
sponding to the connected topology of the second term in the r.h.s. of cq. (45)). 
It is convenient at this point to denote 

G-+{x,y)^{<l>{x) ( t>{ y )) j , (50) 

so that one has 

. = / *Y e i(*i+*s)-" G-+ (y, . (51) 

The definition of the object G^ + (x, y) is identical to the usual definition of the 
— h component of the Schwinger-Keldysh propagators, except, as the notation 
suggests, it is constructed from the classical solution of the equations of motion 
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rather than from the full field operator. Inserting this definition into eq. (48) 
and keeping only the lowest order 16 in the gradients in X, one obtains 

n s (x, P ) = ^lJ-0- 4 G^(x,k)G^(x, P -k) 

+ other contractions . (52) 

Note that in this case there is only one other contraction, that leads to the same 
contribution, thereby transforming the prefactor 1/4 into a 1/2. If the time X° 
at which this is evaluated is large compared to (Q s ) _1 , the classical field that 
enters in the definition of is mostly on-shell, and one can write 

G~+(X, k) w 2ir5(k 2 - to 2 ) f cl {X, k) . (53) 

By analogy with eq. 28, the distribution f c \(X,k) can be interpreted as repre- 
senting the "particle content" of the classical field. As eq. (52) has exactly the 
structure of a 2 — > 1 collision term with on-shell particles of equal mass, it is 
zero because of energy-momentum conservation. 

Therefore, to correctly estimate the magnitude of the source term II s when 
the classical field is weak, one needs to properly account for the slight off- 
shellness of the field Fourier modes. From the equation of motion 



Ql WJ 2 

the off-shellncss of the classical field comes from its self-interactions. The sim- 
plest way to take this off-shellness into account is to use the equation of motion 
in order to write 

**>-^irb/ *•*»-•>■ < 55 » 

and to replace some of the cp's in eq. (48) by the above relation. It is sufficient 
to replace two (fi's in order to lift the kinematical constraints that came from 
the classical field having only nearly on-shell Fourier modes. This substitution 
is straightforward. One obtains, 



U s (X,p) 



\ 2 Ql\ f d 4 k d 4 q 1 



A J J (2tt) 4 (2tt) 4 (fc 2 - to 2 ) 2 



x G cl + (X, q) G cl + (X, k - q) G cl + (X,p - k) + other contractions . 



(56) 



16 At this order, this is equivalent to assuming, from the translational invariance in the 
transverse plane of a large nucleus, that eq. (51) can be replaced by 



(<Hfci)<Kfc 3 )) ~ (27r) 4 <5(fc 1 +k s )G-+(X,k 1 ) 
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This contribution to the source term can be represented diagrammatically 



as 




+ 



(57) 



where the solid lines represent ordinary vacuum propagators (l/(fc 2 — m 2 )) and 
the wavy lines represent the correlation function G~ x + . It is interesting to note 
that this contribution is identical in form to what one would have obtained in 
the collision term of the conventional Boltzmann equation, except that here the 
G h propagators are made up of the classical fields. 

We are now in a position to estimate the power counting of contributions 
to the source term. First, the order of magnitude of the denominators k 2 — m 2 
is Q 2 because the momentum transfer k is of order Q s (and is not particularly 
close to the mass shell). Each G~j + contains a delta function. Two of them can 
be used to perform for free the integrations over the energies k° and q°, while the 
third provides the value of one angular integration variable. We finally obtain 
the estimate 



/* = A~ 2 , n* = QfA -2 and n c \(X) is the spatial density of particles correspond- 
ing to the classical field. The expressions /* and n* correspond respectively to 
the maximal values of f c \ and n c \ can have at early times < Qj 1 )- The argument 
p cannot be specified exactly (in fact, eq. (58) is an oversimplified version of the 
actual formula for II s ), but it is a momentum whose components are of the 
same order of magnitude as those of p, the momentum of the produced particle. 
This is an important point, because as time increases, the support of / c i shrinks 
in the p z direction because of the longitudinal expansion of the system, thus 
making f c \(X,p) decrease as well (while in the center of its support, it would 
stay constant). 

Even if eq. (56) is not valid (say, if the average over j were to generate 
connections among the fields that are not pairwise), the estimate of II s one 
obtains from it has a much wider range of validity. (Eq. (58) is valid even in the 
saturated regime.) We also note that as II s is an inhomogencous term existing 
even when / = 0, its magnitude depends only on the time dependence of the 
classical field 4>{x) through f c \ and n c \. 

5.4 Collision terms 

The estimate of the various contributions to the collision term follow very closely 
that of the source term. Let us start by listing the terms we need to estimate. 
Because of the presence of the background field and of the average over the 




(58) 



where 




(59) 
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external source j, II \.(X,p) can contain topologies that would not exist in the 

vacuum. In fact, II v (X,p) can contain terms that have 0, 1 and 2 loops 17 

before the average over the external source is performed. We will denote by 
Co [/] , C\ [/] and C2 [/] their respective contributions to the collision term. 
Let us start with Co [/] . Diagrammatically, it corresponds to 





Co[f}= — H ¥*—= ■ (60) 



Note that this represents only one of the diagrams that can possibly enter in 
Co [/] . From the experience gained in the estimate of the magnitude of II s (X, p) , 
we can readily see that there must be at least four insertions of the classical 
field for such a contribution to be kincmatically viable when the classical field 
becomes weak and has only near mass-shell Fourier modes. The second equality 
shows one example of the topology one obtains after the average over j. The 
corresponding expression reads 



com 



X 2 Q 2 S \ [ d A k d 4 q 1 



4 J J (2tt) 4 (2tt) 4 (fc 2 - to 2 ) 2 
x G~ + (X, q) G~ l + (X, k - q) G~ l + (X,p- k) + other contractions . 

(61) 

The only difference between this expression and that of II s in eq. (56) is that 
one of the correlators is now replaced by 18 

G-+(X,p) = 2nS(p 2 - m 2 )f(X,p) , (62) 

that involves the distribution f(X,p) rather than the classical distribution 
fd(X,p). From this analogy, we can estimate the magnitude of Co[f] directly 
from that of II s in eq. (58), by substituting one factor / c i or n c \ by respectively 
/ or n. Here n is the spatial density defined from / in the same way as in 
cq. (59). We obtain 

Q 2 S r f n cl (X) \ 2 f(X,p) n cl (X) n(X) U^ph 
Similarly, C\ [/] corresponds to diagrams of the type 





17 Naturally, there are also terms with an even larger number of loops, but these are sup- 
pressed if the particle occupation number is / <C A~ 2 . 

18 This formula for the correlator is only valid as long as the occupation number 

f(X, p) is large compared to one. Its full expression contains 0(p°) + f(X, p). 
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and the corresponding expression reads 



\ 2 Q 2 S \ 2 f d A k d 4 q 1 



(2tt) 4 (2tt) 4 (fc 2 - to 2 ) 2 
x G~ + (X,q) G^ + (X 7 k- q) G~ + (X,p- k) + other contractions . 

(65) 

Here we replace two out of three correlators G~j + by G h ; the power counting 
for this diagram is then 



MjQ n(X)/(A,p) ^ ^n(X) V /ci(*,j>) 
n* n* /* 



/* 



(66) 



Finally, for the 2-loop contribution to the collision term, we have 



C 2 [/] 





(67) 



A 2 Q 2 \ /* d 4 fc A 1 



and 



(2tt) 4 (2tt) 4 (fc 2 - m 2 ) 2 
x G~ + {X, q) G~ + {X, k-q) G~ + (X,p - k) + other contractions , 

(68) 



Ql fn(X)\ 2 f(X,p) 



5.5 Discussion 

Following the power counting in equations (58), (63), (66) and (69), we are now 
in a position to discuss qualitatively the relative magnitude of the various terms 
at different stages of the evolution of the system. An important facet of the 
temporal evolution is that the functions f c \ and n c i are determined once and for 
all from the classical field <p(x) itself. They do not receive any feedback from 
the particle distributions or densities, denoted by / and n respectively, that are 
created in the evolution by the source term II s . The time dependence of n c \{X) 
is driven by the expansion of the system; therefore at times larger than (Q s ) _1 , 
one has 

1^)1 (70) 
n* Q s t V ' 

This reduction of the classical particle density n c \ with time happens because 
the support in momentum space of the corresponding phase-space density f c \ 
shrinks. At a given space-time location X (specified by the space-time rapidity 
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rf), only particles with a matching momentum rapidity y = rj can stay for a long 
time. Therefore, inside its support, f c \ remains constant satisfying 

Ji 1 • ( 71 ) 

Note that at times smaller than (Q s ) _1 , f c \/f* and n c \(X)/n* are also both of 
order 1 because the classical field is completely saturated. 

However, in all the estimates of the previous subsection, f c \ is evaluated at 
some arbitrary location X and momentum p. Therefore, p will eventually fall 
outside of the support of f c \, and f c \ will decrease quickly 19 after that happens. 
For / c i, which comes entirely from the classical field <p, the only time-scale in 
the problem is l/Qs and thus we expect f c \ to start decreasing at times larger 
than l/Qs- 

At early times, r — > 0, the system does not have particles yet and we have 
/ = n = 0. Obviously, in this regime, only the source term II s is important in 
the right hand side of the Boltzmann equation. The corresponding physics is 
that a population of particles, described by the occupation number /, is built 
up from the decay of the classical field. However, these particles are still too 

few to have collisions at a significant rate. Eq. (58) tells us that II s ~ in 
this regime. 

As a rough estimate, if we integrate this source term in the range < t < 
Q^ 1 , we find that the occupation number for particles of momentum p ~ Q s at 
a time r ~ Qj 1 is 

— Y* 1 • {12 > 

At this time, all the components of the momenta of these particles are typically 
of order Q s . Therefore, we also have 

^ = Q ^~1. (73) 

At times around (Q s ) _1 ah the terms in the right side of the Boltzmann equation 
are of equal magnitude. Indeed, in this regime, terms with an arbitrarily large 
number of loops contribute equally to the collision term when / ~ /* . There 
would therefore be an equally large C3[/],C/i[/], etc... In practice, this means 
that one should start using the Boltzmann equation only at later times. 

At later times, r > Q^ 1 , collisions among the particles become important 
and their qualitative effect is to broaden the momentum distribution of the par- 
ticles represented by /, thereby counteracting the effect of the expansion 20 of 
the system. Thanks to these collisions, f(X,p) falls at a lesser rate compared 
to f c i(X,p) (which is not affected by collisions), which eventually leads to the 



19 The precise time dependence of this fall depends on the p z dependence of f c \. To take an 
extreme case, there would be no fall at all if / c i is independent of p z . 

20 In the absence of collisions, / would be affected by the system expansion in a similar way 
to / c i, and its support would shrink like r _1 in the p z direction. 
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dominance of C2 [/] over all the other terms in the right hand side of the Boltz- 
mann equation. When this occurs, our Boltzmann equation is identical to the 
usual one. The detailed mechanisms of this transition between the classical field 
dominated regime and the kinetic regime will be discussed in a future work. In 
particular, it will be interesting to compare, for the QCD case, the temporal 
evolution of the kinetic equation for the glasma with the "bottom up" scenario 
of thermalization [41]. 

6 Summary and Outlook 

In this work, we developed the formalism of Refs. [1,2] for particle production 
in the presence of strong sources to construct a kinetic theory relevant for the 
early "glasma" stage of a heavy ion collision. In particular, we considered for 
simplicity, the dynamics of a <fi 3 theory in the presence of strong sources. Much of 
our discussion however is completely general and could in principle be extended 
to describe the dynamics of gauge fields exploding into the vacuum after a heavy 
ion collision. We showed that the relevant kinetic equation for the particle 
distributions / has the structure of a Boltzmann equation with an additional 
inhomogeneous (/-independent) source term denoting particle creation from the 
decay of the classical field. The collision terms in the Boltzmann equation also 
have novel features. In addition to the usual contribution from the two loop 
self energy, there are 0-loop and 1-loop contributions that affect the particle 
phase space distributions. We outlined the power counting that controls the 
magnitude of the contributions of the source term and the collision terms. The 
temporal evolution of these contributions was discussed only briefly and will be 
discussed in detail elsewhere. 

There are several unresolved issues that should be addressed in future work. 
Primarily, we would like to understand precisely how the derivation here plays 
out in the QCD case. In Refs. [42-45], it was shown that instabilities of the 
Weibel type [46-51] can spoil the bottom up scenario of thermalization. Such 
an instability is also seen in the CGC framework in the explosive growth of 
small fluctuations about the classical background field [52-54] and has a natural 
interpretation as quantum fluctuations about the classical background fields 
on the light cone [55]. A numerical study of instabilities in a field+particle 
framework has been performed [56] but we would like to better understand how 
the effects of such instabilities manifest themselves in the kinetic equation for 
the glasma. It would be especially interesting to uncover whether Kolmogorov 
turbulent spectra [57] arise as a consequence of these instabilities [58,59] and 
whether this phenomenon of "turbulent thermalization" can be accommodated 
in our kinetic framework. 
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A Average over the sources 

We have seen that it is crucial for the validity of the gradient expansion to 
consider quantities averaged over the source j coupled to the fields. We shall 
discuss briefly here how this average can be accounted for in our formalism. 

Let us start from the generating functional for Green's functions of the 
Schwinger-Keldysh formalism 21 Zj[q], for a given configuration j of the ex- 
ternal source. We define it in such a way that the n-point Green's functions is 
obtained by differentiating n times with respect to 77, and then by setting the 
auxiliary source if to zero. From what we have said in section 2, this generating 
functional is related to the sum of all the vacuum-vacuum diagrams by : 

ZjH =e lV s K l3+v\ ) ( 74 ) 

where we have again used a compact notation compared to eq. (9). We do not 
use a boldface letter for the external source j, in order to emphasize the fact 
that it is identical on both branches of the closed time path. 

From this object, it is very easy to construct the generating functional for 
Green's functions that are averaged over some ensemble of external sources, 
with a distribution W[j], as : 

Z[r]} = J [Dj] W[j] e lV sKti+v] ( 75 ) 

In order to see how this average over j can be accounted for in the Feynman 
rules, it is useful to write the generating functional for a fixed j as follows : 

e iV SK li+n] =cxp ^ J d i xV 

x exp (-1 J d*xd 4 y (j(x) + rj(x))G°(x, y)(j(y) + V (y))) , (76) 

where V is the sum of all the interaction terms in the theory under consideration 
(i.e. all the terms of the Lagrangian density that are of degree > 3 in the field). 
In this formula, G°(x,y) denotes the free propagator in the Schwinger-Keldysh 

21 In order to keep the notations compact, we denote by a boldface letter r] the pair r) = 
(v+>V-), where the ± indices refer to the Schwinger-Keldysh closed time path. 
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formalism (as opposed to the full propagator defined in eq. (9)). It is now 
convenient to write the second exponential in the r.h.s. of eq. (76) as the action 
of a translation operator on a functional that does not depend on j, 

exp(-i f d 4 xd 4 y (j(x)+r ) (x))G°(x,y)(j(y)+r,(y))) = 

= cxp^ J d 4 zj(z)-^-^j cx p(~^ J d AxdA V n(x)G°(x,y)r](y) S j . 

(77) 

By inserting this formula in eq. (76), and then in eq. (75), we obtain the following 
expression : 

Z H = {J M W W e i/cd4zj(z) ^} 

x exp (i J^d 4 xV (j^jj) exp (-^ J^d i xd i y r 1 (x)G a (x,y)r 1 (y)^ . 

(78) 

The terms on the second line are nothing but the generating functional for the 
same theory without any external source (since it does not depend on j). As 
we can see, the effect of the average over the external source j is to bring a 
prefactor which is a certain functional of the operator S/drj. Such a term can 
be interpreted as additional couplings among the fields, since one can always 
write : 

[Dj] W\j] e*/c-*»JW4iy} ee exp (i J^x U (^y)) • (79) 

What this derivation makes obvious is that, for calculating averaged quantities 
over the ensemble of external sources j, one can forget the external sources 
altogether, and include additional vertices to the theory 22 , as prescribed by 
eq. (79). Note that this is equivalent to calculating a quantity in an arbitrary 
j, and then reconnecting all the j's among themselves in all the possible ways 
permitted by h^W-^j]). 
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